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SUmRY 


An analytical technique for the prediction of fan blade flutter is 
evaluated by utilizing firstostage fan flutter data available from tests 
conducted at the National Aeronautics and Space Administration (NASA) » Lewis 
Research Center on an advanced high performance engine. The formulation 
Includes both aerodynamic and mechnlcal coupling among all the blades of the 
assembly. Mlbtunlng is accounted for in the analysis so that individual 
blade Inertiasi frequencies, or damping can be considered. Airfoil stability 
is predicted by calculating a flutter determinant, the eigenvalues of which 
indicate the extent of susceptibility to flutter. 

When b lade- to-b lade differences in frequencies are considered in the 
analysis a stable system is predicted for the test points examined in this 
report. For a tuned system, however, it was found that torsional flutter can 
be predicted at a limited number of interblade phase angles. An examination 
of these phase angles indicated that they were "close" to the condition of 
acoustic resonance. For the range of Mach numbers and reduced frequencies 
considered here, the so-called subcrltlcal flutter cannot be predicted. 

The essential influence of mechanical coupling among the blades is to 
change the frequencies of the system with little or no change in damping. 

But, aerodynamic coupling together with mechanical coupling could change not 
only frequencies but also damping in the system, with a trend toward 
instability. 
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INTRODUCTION 


An Important factor which influences the development program of any 
advanced engine is the susceptibility of rotor blades to flutter Instabilities. 
In order to understand more fully the mechanism of fan blade flutter Instabll-* 
Ities and the parameters that influence flutter) several tests on an advanced 
engine were conducted at NASA Lewis. These tests covered a wide range of 
realistic operating conditions appropriate to a modern high performance jet 
engine. Considerable data were generated for the component under study (i.e.t 
the part'-span-’Shroud fan assembly). Documentation of the steady**state aero- 
dynamic conditions prevailing at the several test points, compilation, and 
analysis of flutter data may be found in reference 1. Detailed analysis of 
aerodynamic data obtained from high response pressure probes and optical 
probes are presented in reference 2. 

The nature of flutter observed in all these tests was subsonic stall 
flutter and a study of the flutter mode of vibration revealed the presence of 
several harmonics in a nearly pure torsional mode. Frequency measurements 
made on each blade fully restrained at its shroud revealed blade-to-blade 
differences. Considerations of all these aspects of the problem pointed to 
the need to develop an analytical model which can account for blade mistuning 
in the analysis. This report summarizes the results of a preliminary effort 
to model the flutter characteristics of a mistuned rotor. 

In as much as the observed flutter mode was concluded to be essentially 
an above shroud first torsional mode, the analytical model represents each 
blade as a torsional oscillator. The oscillators are assumed to be mounted on 
a ring representing the shroud. The formulation allows for blade-to-blade 
differences in any or all of the following; frequency, inertia, damping. The 
unsteady aerodynamic coefficients used in this study are based on the analysis 
due to Smith (reference 3). The problem is formulated along the lines out- 
lined by Whitehead (reference 4) and leads to a flutter determinant, the 
solution of which yields the eigenvalues at which flutter is likely to occur. 
The real parts of the eigenvalues represent the frequencies at which flutter 
may occur and the corresponding imaginary parts represent the level of 
damping. 

In this report, the technical approach, along with a discussion of its 
principal features, is outlined first. The application of the analytical 
model to four test points is discussed next, followed by a discussion and 
summary of the results covering the conclusions drawn from the comparisons 
made between the analytical model and test results. The detailed derivations 
are shoxm separately in the APPENDICES. 
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NOMENCLATURE 


A 



Cm 

•(c) 


c 

E 

(E) 

f 

I 

K 

M 

(M) 

n 

N 

P 

S 

t 

u 

U 


Aeroelastlc moment 

Matrix of aerodynamic effects 

Unsteady aerodynamic moment coefficient 

Matrix of Cm 

Blade chord 

Modulus of elasticity 

Matrix of Interblade phase angles, see Eq. (11) of APPENDIX A 

Force between blade and support structure 

Mass moment of Inertia of a blade/unit length or 
moment of Inertia of ring cross section 

Matrix of blade Inertia ratios 

Blade stiffness 

Mach number 

Matrix of mechanical effects 
Nodal diameter (hamonlc number) 


Number of blades 
Frequency parameter 



Gap between blades 



w 


Aerodynamic modal twist 

Steady fluid velocity relative to blade 

Amplitude of vibratory rotation of support structure 
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y Amplitude of vibratory rotation of support s|;rdcture 

2 y-w 


a 

/3 

y 

8 

h 

V 

V 
X 

e 

p 

c 

Subcripts 


Dynamic Influence coefficient 
Interblade phase angle 

Individual blade frequency departure parameter 
Logarithmic decrement 

Ratio of aerodynamic inertia to blade inertia 

Frequency parameter 1 — — 

Elastic axis position expressed as a percentage of chord from the 
leading edge 

Circular frequency of vibration 
Reduced frequency ctu /U 

Blade stagger relative to tangential plane 
Air density 
Damping ratio 



AR Acoustic resonance 

a Twist 

F Force 

I Imaginary part 

M Moment 

o Leading edge position (or reference blade) 

*1 Translation 

R Real part 

T Tangential 
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TECHNICAL APPROACH 


Overview 

The approach adopted here considers the differential equations of motion 
governing the vibration of a cascade of blades in a uniform flow field. Thust 
the forces acting on each member of the cascade include aerodynamic forces in 
addition to mechanical (inertial, restoring, and damping) forces. The aero- 
dynamic forces acting on the system not only influence the vibratory motion of 
the blades but are, in turn, influenced by vibration. These aeroelastlc 
forces are generally dependent on the prevailing steady-state aerodynamic 
conditions, cascade geometry, and the mode of vibration. For a complete 
description of the aerodynamic force it is also necessary to know the location 
of the elastic axis. As each blade is a member of an assembly of blades 
vibrating in a flow field, the vibration of each member influences every other 
member and thus all blades are aerodynamically coupled. In addition, the 
blades can also be coupled mechanically through the disk on which they are 
mounted or through the shrouds at v;hich each blade comes into contact with its 
neighbors due to blade untwist at speed. 

The above considerations lead to a formulation which describes the motion 
of a system of blades coupled aerodynamically and mechanically. By virtue of 
mistuning in which one or more characteristic (frequency, damping, etc.) of 
each blade is "slightly" different from those of a reference blade, the system 
of equations leads to a flutter determinant, the solution of which yields the 
possible frequencies of vibration at which the system is susceptible to 
flutter. 

The construction of the mathematical model proceeded on the basis of 
representing only the above-shroud-torsion vibration mode of an advanced fan. 
Thus the analysis is restricted to modeling the blades as linear torsional 
oscillators supported on a flexible "shroud ring". The validity of the model 
and the associated computer program was established by making comparisons with 
appropriate published analytical results. Upon validating the model, several 
flutter test points were chosen from available data and analyzed to determine 
if the analytical results indicated trends observed in tests. 
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Analytical Model Development 

The details of the analysis are presented in APPENDIX A. The complex 
eigenvalues are sought as a solution to the matrix equation 

( <» 

where the matrix Q ^ - I represents mechanical i=ffects and the matrix 
(Q) (E) (C) (E)“l represents the aeromechanical effects; (O) Is a matrix 
whose elements Include the dynamic Influence coefficients of the support 
structure; fx represents the Inertia ratio; and X is the reduced frequency. In 
terms of Whitehead's notation the eigenvalues p are calculated from using 


the simple relationship 



Unless otherwise indicated, the results presented in this report show the 
values of p rather than l/. The real part of p represents the departure of the 


( 2 2 
D = 

and the Imaginary part of p represents damping in the system in terms of the 
loss factor. The latter is a commonly used measure of structural damping and 
is defined as the ratio of specific damping energy to’ the maximum strain 
energy per unit volume of the vibrating body. 

Some of the principal assumptions inherent in the analysis are worth 
noting. The dynamic influence coefficients of the "shroud ring" which are 
derived in APPENDIX B are dependent on th(, frequency of vibration of the 
entire system. Thus th «2 ring is characterized at a reference frequency in much 
the same way the disk was characterized in reference 5- Further, the manner 
in which X enters into the analysis is such that the calculations are made for 
a fixed X . 

The analysis and the associated computer program are developed in such a 
way that the unsteady aerodynamic coefficients needed in the calculations may 
be obtained from any available analysis and used as input. However, the 
results reported .here, are all.bssed on the use of the complex coefficients 
calculated by the method given by Smith (reference 3). Smith's analysis is 
based on the assumption that the blades can be represented by flat plates at 
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aero mean Incidence » that the flow le subsonic » isentroplci and two-dimension- 
al and that the Kutta-jpukowski condition at the trailing edge is satisfied. 


Calculation Procedure 

The unsteady aerodynamic coefficients were first calculated for a given 
cascade geometry and these were transferred to the appropriate elastic axis 
position. The extent of mlstunlng in terms of departure of blade frequency 
and/or damping and/or blade Inertia from those of the corresponding mean 
values are used as input. To Include mechanical coupling among the blades 
through 'the supporting "ring", the dynamic influence coefficients appropriate 
at the reference frequency were calculated and used in computing the matrices 
/fand At Eigenvalues and eigenvectors! were then computed using the method 
of approach available in the EISPAK package (see reference 6). 

For a single degree of freedom representation, a positive sign on the 
imaginary part of the aerodynamic coefficient Implies negative aerodynamic 
damping and, therefore, for a tuned rotor no further calculations will be 
needed to determine which of the test points were likely to indicate flutter. 
The condition for incipient flutter is therefore obtained by setting the 
imaginary part of the aerodynamic coefficient to zero, i.e., 

A solution of this quadratic in 7j provides the positions of the elastic axis 
appropriate for flutter to occur for prescribed values ofXandjS. However, 
when mistuning is included, a complete solution of the flutter determinant is 
necessary so that only the nature of the complex eigenvalues for the system 
can determine the extent of flutter. 

Analytical Model Verification 

A series of results was generated for a twelve-bladed system in order to 
validate the accuracy of the computer program and to illustrate the various 
features of an aeromechanically coupled assembly of vibrating blades. 

Figure 1 illustrates the eigenvalues for the undamped and tuned blades 
supported on a flexible ring. Only mechanical coupling effects are shown in 
figure 1. As the support structure does not distinguish between forward and 
backward traveling waves, only six (N/2) distinct eigenvalues are calculated. 
The ring was characterized at a reference frequency (Fq = 1098 cps) . As 
expected, some of the eigenvalues for the mechanically coupled system are 
higher than those of the reference frequency and some lower. The extent of 
departure depends entirely on the proximity of the ring eigenvalues to the 
frequency at which the ring is characterized. 



On the other handt If the same system is coupled through aerodynamics » 


li the eigenvalues are unsymmetrlcally distorted as shown In figure 2 and result 

|i In a locus, the real part of which represents the frequencies of vibration of i 

I; the system and the Imaginary part represents the damping. The asymmetry 

i; (l.e., the frequencies of the forward and backward traveling waves being some- 

k what different from each other) Is due to the dependence of the aerodynamic , 

|i coupling forces on the direction of the wave (l.e., the sign of the Interblade 

ii phase angle). > 

I The results of mlstunlng alternate blades Is shown in figure 3. Figure 3 

also Illustrates the manner in which the locus of the eigenvalues shifts with 
(a) the addition of mechanical damping, and (b) the removal of all aerodyiiamlc 
coupling. From figure 3 it is clear that the locus of eigenvalues shifted 
upward (toward more stable condition) Indicating a correct trend. Further, 
with the removal of all aerodynamic coupling, the locus reduces to clusters 
of points. The magnitude of the imaginary part of the eigenvalue (l.e., 

0.00A4) is 8 /tt so that 8= 0.0137 which agrees with the logarithmic | 

j: decrement input. j 

The cascade parameters used to generate the loci of figure 2 are the ♦ ’ 

same as those of reference 4. The unsteady aerodynamic coefficients were 
calculated for interblade phase angles )S=2Tm/N from n=?l to N. The departure 
i of each blade frequency from the mean ( y ) and the damping distribution (8) 

were obtained by enlarging figure 18 of reference 4 and picking out the 
numbers from the enlarged figure (sec TABLE I). These were then used in cal- 
culating the effects of mistuning and the results are shown in figure 2. 

The eigenvector corresponding to the least stable eigenvalue is shown in j 

figure 4. The lack of corrcjlation at one point in figure 2 is attributed to I 

the fact that /3=0 ( or 2ir) is specifically excluded in the computer program 1 

used to calculate the unsteady coefficients. The combination M=0,jS =0 is of 

little Interest in practice and therefore no further studies were made. The 

strength of each harmonic present in a mistuned mode is calculated as a part i 

of the computer program as shown in figure 5. As observed by Whitehead in 

reference 4, the eigenvector has a predominantly two nodal diameter pattern, 

as can be seen from figure 5. 

Figure 6 shows the effects of both aerodynamic and mechanical coupling on | 

1: the locus of eigenvalues. Due to the addition of mechanical coupling, the | 

locus of eigenvalues does not shift uniformly and the frequencies correspond- . 

! Ing to an interblade phase angle may increase, decrease, or experience little 

k or no change at all. Comparison of figures 6 and 1 indicates the nature of 

f frequency shifts shown in figure 6 corresponds closely to those shown ^ * 

figure 1. Thus, the eigenvalues corresponding to nodal dimae ter patterns -3 

and -4 in figure 6 are affected most and in a manner similar to that shown in 
figure 1. The mechanical coupling corresponding to figure 6 was obtained for a 
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stiff and "oversized shroud ring"* A more realistic case corresponding to a 
more flexible coupling among the blades will be considered later In connection 
with the analysis of a 38 bladed rotor. 

Reference to figure 7 confirms the nature of the eigenvectors typical of 
a mlstuned system. The wavlness present In the phase angle confirms the non- 
constant nature of the Interblade phase angle; an Important feature of a mls- 
tuned system. Further, an eigenvector corresponding to an eigenvalue for a 
mlstuned systetn contains In It contributions from several modes of the tuned 
system. 


Mechanical Coupling Among Blades 

Sines the vibration modeling is restricted tp considerations of above- 
shroud torsional mode only, the shroud is assumed to serve as a support for 
the blades. Thus the deformation of the shroud couples all the blades. The 
shroud is modeled as a coutlnuous ring whose dynamic influence efficients were 
calculated at a frequency of 1098 Hz; this frequency being the average of the 
38 blade frequencies. The latter were obtained from reference 1. The 
influence coefficients exhibited the expected characteristics which included 
(a) identical numbers on the principal diagonal, (b) complete symmetry, and 
(c) identical numbers on each subdiagonal. For a 38 bladed rotor, calculations 
were made using the actual dimensions of the "shroud ring" supported by 
springs in a plane perpendicular to that of the ring. The stiffness of the 
springs (139.931 N/m/radian) was obtained from an available NASTRAN run for 
the zeroth harmonic of the fan assembly. 


APPLICATION OF ANALYTICAL MODEL 


The test points at which flutter was observed in the advanced fan are 
marked on the fan map shown in figure 8. The analysis procedure and the 
associated computer program were applied to calculate the susceptibility to 
flutter for four test points (i.e., Nos. 118, 123, 125, and 127). The 
steady-state aerodynamic conditions for these test points are fully documented 
in reference 1. The first step was to calculate the unsteady aerodynamic 
coefficients for each of these points for a representative section of the 
blade and examine their behavior. TJie tip section was chosen as the 
representative section in this study. 

Figure 9 shows the loci of the aerodynamic coefficients. An examination 
of the coefficients for points 118 and 125 shows the extent of sensitivity of 
the unsteady aerodynamic coefficients to the parameters X and M. Although the 
difference in X and M between these test points is only 0.01 and 0.002, 
respectively, test point 125 indicates flutter (Cj >0) whereas test point 118 
does not. While computing these coefficients it was noted that their magni- 
tude varied somewhat with the number of collocation points used on the airfoil. 
An initial study of this was made for test point 125 using four collocation 
points (Instead of five used in all other calculations) and the results are 
shown in figure 10. It is noted that even for four collocation points, the 
results are within experimental accuracy. 

At the Mach numbers (MaO.8) and reduced frequencies (Xa l.8), it was 
noted that the aerodynamic coefficients formed a nearly continuous curve 
except for a few (three in the cases of test points 118 and 125, for example) 
coefficients which appeared separate and distinct from the rest. The inter- 
blade phase angles corresponding to these particular coefficients were checked 
to see if they are "close" to either /8®satisfying the acoustic resonance 
condition. The latter is discussed in some detail in reference 7 and the 
interblade phase angles can be calculated from the "cut-off" condition to be 

^ar = -f ( ^ 6 ± - M^cos^e ^ 

The regions of flutter in the (17, /3) domain are obtained by solving Eq. (4) and 
are shown in figure 11 for certain cascade conditions which Include those of 
test point 125. Values of /S calculated from Eq. (4) are also shown in figure 
11. It is clear that under the assumptions discussed earlier, the nature of 
flutter implied for the test conditions of points 118 and 125 is acoustic 
resonance. As can be seen from figure 11, this phenomenon occurs over an 
extremely narrow region of interblade phase angles corresponding to backward 
tfaveling waves . The physical significance of this phenomenon in the case of 
a real machine is not clear. 



I Individual blade frequencies were available for the fan assembly under 

i Study and these had been measured (reference 1) with each blade fully 

[ I restrained at the shroud. The nature of variation of these frequencies Is 

i I * shown In figure 12. Using these frequencies, the mlstuned system was 

i i analyzed first with aerodynamic coupling only and the results are shown In 

! f' figure 13. VIhlle these results show the system to have positive aerodynamic 

I I * damping, It Is essential to bear In mind that (1) the aerodynamic coefficients 

f I are very sensitive to changes In X and M, as can be seen from a comparison 

I of the coefficients for test points 118 and 125 (see figure 9), and (2) a 

I mlstuned mode consists of many Interblade phase angles. Thus It Is likely 

I that with slight changes In the magnitudes of Important parameters, some of 

: I the Interblade phase angles coincide with those calculated earlier In the 

I acoustic resonance condition. 

II 

f i For the eigenvalue corresponding to the lowest aerodynamic damping In 

I figure 13, the mode shape was calculated and Is shown In figure 14. The 

\ features of this mode are similar to those of the twelve-bladed system shown 

I in figure 7 and include nonconstant Interblade phase angle around the rotor. 

' \ The mixture of "tuned modes" in the eigenvector calculated at the least 

positive eigenvalue is clear from figure 15. These calculations show pre- 
I ’ dominant backward traveling waves of nodal diameter patterns 5 and 6 in con- 

trast to a mixture of several forward and backward harmonics Kurkov ^..reference 
j 8) obtained via analysis of optical displacement data recorded during flutter 

(see figure 16). 

Figure 17 shows the effect of including aerodynamic and mechanical 
coupling on the eigenvalues of the system. In obtaining these results the 
"shroud ring" was modeled by using the actual dimensions of the shroud. Two 
features of these results are worth noting: (1) the lowest aerodynamic damping 

level dropped from 0.0196 (see figure 13) to 0.0056 upon Including mechanical 
coupling; and (2) the frequency at which damping Is the lowest changed from 
1175 Hz (7 percent higher than the average blade frequency) to 994 Hz (9.5 
percent lower than the average blade frequency). 

The mode shape corresponding to the lowest value of aerodynamic damping 
in figure 17 is shown in figure 18. The strength of several harmonics present 
in the mode is shown in figure 19. Clearly, a five nodal diameter backward 
traveling wave is the predominant component. 

It may be noted that the position of the elastic axis (77= 0.44) was 
chosen on the basis of optical displacement data obtained during flutter. 
Reference 8 discusses the details of the method used. In order to determine 
the sensitivity of the aerodynamic coefficients to small changes in 17, it was 
decided to calculate the coefficients for M = 0.795, X= 1.79, and 17 = 0.40 
i and compare these with the locus obtained for conditions prescribed for test 

j point 125 (see figure 20). 
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DISCUSSION OF RESULTS 


Among all the factors that enter Into the calculation of aerodynamic 
damping, the most significant is the aerodynamic theory from which the unsteady 
forces acting on the blade are calculated. The theory chosen for use In this 
study Is based on assumptions of a flat plate cascade and no loading effects 
are included. 

It would be appropriate to repeat these calculations with alternate 
aerodynamic coefficients. For a given cascade geometry the calculations were 
found to be very sensitive to changes In Mach number and reduced frequency 
but less so to Inertia ratio. The reduced frequency and Mach number vary 
from root to tip but considerations of a single representative section such 
as the tip section does not take Into account these variations. Needless to 
say, a realistic solution should consider modeling the entire blade from root 
to tip. However, with mlstuning Included, the problem size can be large 
unless aerodynamic modal representations are used. In much the same way 
structural dynamic modal representations are made by use of generalized mass/ 
inertia and stiffness. 

The importance of the phenomenon of acoustic resonance arises from the 
Implication that at the interblade phase angles satisfying Eq . (4) , the waves 
carry energy in a purely tangential direction. The prediction of this type of 
flutter for the fan considered here Is similar to that made in reference 7. 

The range of Interblade phase angles at which acoustic resonance flutter can 
occur is very narrow and may not correspond to any of the discrete phase angles 
calculated for a finite number of blades In a cascade. However, as the 
unsteady pressures acting on blades Include a variety of harmonics which are 
not necessarily Integral multiples of the first harmonic, only a detailed 
analysis of pressure data can answer the question of the relevance of acoustic 
resonance phenomenon to fan blade flutter. 

Mechanical coupling among the blades could change both the real and 
Imaginary parts of the eigenvalues, The extent of such change appears to 
depend upon the level of stiffness of the support as can be seen from a 
comparison of figures 6, 13, and 17, The dimensions of the ring used in 
obtaining the results shown in figure 17 lead to a much more flexible ring 
than the one used in obtaining the results shown in figure 6. The model used 
in this report characterizes the "shroud ring" at any one frequency and this 
was chosen to be the average of all blade frequencies. A drop in frequency 
of nearly 10 percent from the average suggests a support somewhat more 
flexible than is present in the real structure as the observed flutter 
frequency was in the vicinity of 1078 Hz, a mere 2 percent below the average. 
While the size of the ring cross section could be adjusted to obtain a 


smaller drop In frequency, it was thought that other assumptions in the 
analysis should be studied more thoroughly before any revision of the shroud 
model Is considered. 

The analytical model developed here considers only one degree of freedom 
for the blade motion and thus ignores any coupling that tray be present among 
the vibratory modes. Recently Bendlcksen and Friedman (reference 9) reported 
that the Influence of coupling on blade flutter calculations Is significant. 
For assemblies which have closely spaced modes. It appears that the inclusion 
of coupling among the modes of a blade Is an appropriate next step In the 
development of a flutter prediction system. 


SUMMARY OF RESULTS 


For the cascade conditions appropriate to the test points studied In 
this report, the aerodynamic theory that was used cannot predict subcrltlcal 
flutter. Under the assumptions of a tuned assembly, the imaginary part of the 
aerodynamic coefficients does indicate flutter for a limited number of inter- 
blade phase angles, but these Interblade phase angles are "close" to those at 
which the so-called acoustic resonance is predicted. Upon using the individual 
blade frequencies and solving the mlstuned system with aerodynamic coupling 
only, the results show a stable system. Eigenvectors calculated for the 
mlstuned system showed the presence of several harmonics in each mistuned 
mode. Inclusion of both mechanical and aerodynamic coupling, in the solution ' 
of the eigen problem, influences not only the frequencies but also damping in 
the system with a trend toward instability. 
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APPENDIX A: FOBMUUTION OF ANALYTICAL MODEL 


Governing Equations of Motion 
The vibration of blade Is given by 


( 1 ) 


where 


2k= Yk *Wk 


( 2 ) 


Let 

fk'CkZk+KkZk (3) 

where Mass moment of Inertia of blade blade/unit length 
Kjj = k*-^ blade stiffness 

C{( » Mechanical damping associated with k^^ blade 
yjj = amplitude of vibratory motion of k*"^ blade 

t*Vk 

wk = amplitude of vibratory motion of support structure at k' blade 

Ak = aeroelastlc forces acting on k*"" blade. 

The vibratory motion around the support structure (shroud "ring" or disk) at 
any blade location may be written as 

N 

WK=Iak4®fi (4) 

f-i 

Where Qt\{ represents the vibratory motion of the support structure at k due 
to a unit stimulus at i. 

Clearly, , where w is the frequency (of the forcing function) 

and N Is the total number of blades. 


Vk “ 


( 6 ) 


and with y: y « yk^w^- Cthk^h where the superscript b refers to blade 
Thus, 

Xk' I®kf^« (7) 

where * . . h 

«krOki + «kk 

Equation (1) may then be rewritten as 

-w2(^)(a)jf|-W+lf|=0 (8) 

where ( ^ ) Is the Inertia matrix 

(cd is the matrix of dynamic influence coefficients of the support 
structure and blade 

|a[ is a column matrix representing the aeroelastic moments. 

Aeroelastlc Moments 

Let Uj. represent y^ in the rth mode. The total amplitude of the blade 
is the sum of all complex amplitudes from all contributing modes. 



Where jSf is the Interblade phase angle corresponding to the r*-^ mode 
(i.e., and k/Sr =r(3[^)• Thus, 

|y|=(E)|u| (10) 


where (E) represents the matrix consisting of interblade pahse angles, i.e.. 


(ID 



The aeroelastlc moment acting on the bladei due to motion In the mode 
alone can be written 4s 

= 7r/ouVcM,Ure“‘^^e'‘^’ (12) 

where p Is the air density 
U Is the air velocity 
c is blade chord 

CHj. Is the unsteady aerodynamic moment coefficient corresponding to a 
cascade of blades oscillating In the r^^ mode. 

Let 17 represent the elastic axis position expressed as a percentage of chord 
measured from the leading edge of the airfoil. Then, 

((-MO - (Cmo)q “ 77 (^Fq)o " (Cfq)q (13) 

where Cma represents the moment coefficient with 77 » 0. 

Cpa , Cpq represents the unsteady aerodynamic force coefficients due to 
vibratory twist a and displacement q, respectively. 

The total contribution to the aeroelastlc moment acting on the blade, from 
all modes can be written as 

|A| = 7r;oU^V‘-*(E)(C)(E-')(a)|f| (14) 

Therefore, the governing equations of motion In terms of the force jfj 
transmitted to the supporting structure may be written as 

[-a.2(;)(a) + 7T/DU^C^E)(C)(E)"'(a) + l]|f( = (o| ( 15 ) 
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Th« element of ( a ) may be shown to be 

Co} 


(16) 


W ^ ^ to) 

where 9o and Wp refer to the Inertia and frequency of the reference blade. 
The element of (a) is simply a^j 


The governing Eq. (15) can then be vewrltten as 

4'' 2 

)(£■')(«) + 1 


^ f ) 


|f| ' M 


(17) 


The elements ofoare as shown in Eq. (16) and (^) “ ^0 ’ 

-^(EUcUE)''(a) +l]{tl = |o| 

f, 2 . 

The elements of 6 are those shown in Eq. (16) multiplied by Jq(jj and fx- ^ 
representing the ratio of aerodynamic inertia to blade inertia. 

For the sake of clarity, the elements of a are shown below. 

^H*2i 5o ui^ 

Okj = .^0 V Okj® 

^r» t^r» represent the ratio of blade inertia and frequency etc. (i.e., 

Inertia of blade/ inertia of reference blade, frequency of blade/f requency of 
reference blade, etc.). Multiplying throughout by u;o^/<ir, Eq. (19) can be 
written as 


fllkk = 


O)* 

(ti. 


(19) 


{- 


9)A) - -^(E)(c)(E)''(a)+ ^(l) {fhlol 


( 20 ) 


Clearly, the elements of a are the same as those shown in Eq, (19) except for 

the quantity w which Is divided out. Let Wq24i|2 b (i+i/) and Let^» 

(S)(^) " I represent mechanical coupling effects, and (/I) ■= //./X^i(E) (C) <E”^) 
represent aeromechanical coupling effect. Then Eq. (20) can be shown to 
reduce to 



Equation (21) represents the required eigenvalue formulation. 

Clearly, the solution v provides the frequencies (real part) and damping 
(Imaginary) at these frequencies. If p Is defined In accordance with 
Whitehead's notation, as 




0 ) 2 


o 


then p 
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APPENDIX'S: DYNAMIC INFLUENCE COEFFICIENTS FOR "SHROUD RING" 


The shrouds on the blsdes are assumed to £orm a continuous ring . The 
derivation of the dynamic Influence coefficients Is based on the following 
assumptions: 


1. The ring Is circular, homogeneous, and Isotropic. 

2. The ring Is supported on an elastic foundation. 

3. Only out-of-plane small (linear) vibrations are considered. 

4. Effect of rotary Inertia Is Included. 

5. Effect of shear deformation Is neglected. 


The shroud ring model and the forces and moments acting on an element of 
the ring are shown in figure 22. Referring to figure 22, and assuming that 
the elemental angle d0 Is small, the force and moment equilibrium equations 
can be written as 


js" = kw - f + m — =- 

de 


(la) 


dM^ 

de 


+ Mb = 0 


(lb) 


dMe 

dd 


mr* d^ 



(ic) 


Here, w(fl, t) Is the deflection, V Is the shear force, M,j, Is the 
torsional moment. Mg Is the bending moment, m Is the mass per radian, k Is 
the stiffness per radian of the support, R Is the radius of the ring, r Is the 
radius of gyration of the cross section of the ring, f Is the external force 
per radian, and F Is the external moment per radian. The governing equation 
of motion Is obtained by differentiating Eq. (Ic) with respect to 0, using 
Eqs. (la) and (lb), and noting that 


Mb 


El d^w 

302 
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The equation of motion is 


El . a*w mr* a*w .„,4.iL (2) 

•rS- ' 

The appropriate boundary conditions are: 

d^w 

1. w(o,t) =w(27T,t)i — ( 0 ,t) = — (27T,t) (3a 

2. For a single external force. f|;d6. acting alone atd^9\^t 

|^(9..n = o (3 e 

3. For a single external moment, Fj^dd, acting alone at 3=0 

w(3,^,t) =0 

Free Vibrations 

The problem of free vibrations is obtained by setting fsO and FsO in 
Eq. (2). The solution which satisfies the boundary condition, Eq. (3a), can 
be written as 

w(3,t)= 


where 


i>n = On COS n3 


mR(R2+ r^) 


[Eln^n^-D + kR®] , n = 0,l,2,..., 


and 0=3 is arbitrary. The coefficients, a^^, in Eq. (5) are chosen such that 


= MS ^ (r*-i 
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I 

■T 


Oq® (27rmR) 

■ 

0„= [!^(R^n'r')] 




t for n"0, 


, for 02^1 • 


(7a) 


(7b) 


Forced Vibrations 

In view of the boundary conditional Eos. (3b) and (3c), the solution of 
the forced vibration problem is written as 

w = Wf+Wp (8) 

where Wf is the solution for the case of no external moments (l.e., FsO) and 
wp is the solution for the case of no external forces (l.e., f = 0). 

For the case of FsO, the solution can be written as 

w, = 5 , (9) 

where, 

= QnCOS n(a-®K). (10) 


It is noted that the assumed solution satisfies the boundary condition, 

Eq. (3b). 

By substituting Eq. (9) into Eq. (2), multiplying by and integrating 

from 0=0 to 6= 2 tt , one obtains 

i?n + ^n^‘)7n= R/ f . (11) 

,0 

The external excitation is assumed to be harmonic, so that, 

f = f e'"* , (12) 


where lu is the excitation frequency, and f is, in general, complex. 


Further, assuming that u> Is not equal to Xn» the solution to Eq. (11) Is 
written as 


By using Eqs. (9). (11). (12), and (13), one obtains 


where , 




iuA 


Cn= (Xn“W )"' 


(13) 


(14) 


(15) 


Following a procedure similar to the one described above, the solution 
for the case of fsO can be obtained as 


Wp 


where, 




= -On sin r\{6^~6^^) , 


.iwt 


(16) 


(17) 


and F which is, in general, complex is defined by 

A 


F= F e 


iuit 


(18) 


It is noted that Eq. (17) satisfies the boundary condition, Eq. (3c). 

The solution of the forced vibration problem, Eq. (8), can be written as 


w(0,t) =w(e)e‘‘^’ 


(19) 


where. 


W = R I + Z • (20) 


Differentiation of Eq. (19) leads to an expression for the slope, defined by 


where* 


0(8,0 =0(8)e'"' 


( 22 ) 


and 


0 “ = I « « +Tr? tI°- -Mr <1® 


(23) 


Discretization 


To obtain the dynamic influence coefficients for use in the matrix 
formulation of the mlstunlng problem* the solutions for the deflection and 
slope of the shroud ring have to be discretized* 

Let* N be the number of blades* k, m* or j the blade number* n the mode 
number, 0»O be fixed at blade number k»l, andAfl= (27r/N). The discretized 
version of Eq. (20) can be written as 




(24) 


where, 


Cnn ” Cn » 


<^Sn = On COS n (0„ -0h) . 


(25a) 


^mn = - On Sin n(0m-flk). 


(25b) 


f oiT ii“0 p X f 2 • • • Slid Tn j lc"X ^ 2 ^ 3 * • • 


Letting, 


(-) = 

\de /j 


2A0 'J 


• A A . 

(F: + ,-F,.|) . 


Om " ^ ® > ®m “ ® 


in Eq. (24), and after some algebra, the matrix version of Eq. (24) can be 
written as 
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|w[ = [a‘](q| + [q 2] iQf 


(26) 


A 12 

where w, q» and Q are N<*viictors» and a tOe are NxN matrices, whose elements 
are given by 


®mk ' Z Cnn On 

nsQ 


(27a) 


QmK " ’^mn Cnn ,Z ^ • 


(27b) 


m refers to the row, k refers to the column, and 

0 10 0 


[E] = 


-1 


-1 


0 -1 


0 


-1 0 


-1 0 


(28) 


L 10 0-1 0 J 

Following a procedure similar to the one described above, the matrix 
version of Eq. (23) can be written as, 


\$\ » [a’]{q( ♦[a-ltol 


where. 


®mk ■ X ^ tmn ^nn On 
nsO 


(29) 


(30a) 


O mk = - ^Z^ 0 mn Cnn X^ ^ jk • 


(30b) 


Equations (26) and (29) can be combined to give the dynamic Influence 
coefficient matrix as, 

\M HJ 
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TABLE I 

DEPARTURES OF INDIVIDUAL BLADE FREQUENCIES (/) AND 
DISTRIBUTION OF DAMPING ( 8 ) USED IN MODEL VERIFICATION 
Reference 4, figure 18 


k 


S^/tt 

1 

0 

0.0025 

2 

-0.0055 

0.0040 

3 

-0.0015 

0.0020 

4 

-0.0008 

0.0026 

5 

0 

0.0025 

6 

0.0013 

0.0022 

7 

0.0010 

0.0029 

8 

0.0002 

0.0052 

9 

0.0002 

0.0052 

10 

0 

0,0056 

11 

0.0019 

0.0055 

12 

0.0034 

0.0055 
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Figure 2 Effect of Mi! 
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Figure 3 Effects of Mistuning on Blade Flutter 


(SEE FIG. 2 FOR THE LOCUS OF EIGENVALUES) 
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Figure 4 Eigenvector Associated with Eigenvalue p=- 0.0872 — i 0.0123 
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Figure 6 Eigenvalues for a 12 Bladed Tuned System 
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Figure 8 Fan First Rotor Subsonic Stall Flutter Boundaries 


TEST POINT 118 X=1.78 M = 0.797 TEST POINT 125 X=1.7S M = 0.795 



Figure 9 Unsteady Aerodynamic Coefficients 
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Figure 10 Unsteady Aerodynamic Coefficients (Four Collocation Points) 
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Figurell Acoustic Resonance Flutter 
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Figure 12 Departure of Individual Blade Frequencies from the Mean 
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Figure 17 Mistuned Assembly with Aerodynamic and Mechanical Coupling (Test Point 125) 
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Figure 19 Strength of Harmonics Present in the Mistuned Mode at p = —.1804+1.0057 



Figure 20 Unsteady Aerodynamic Coefficients 





